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The  objective  of  this  program  was  to  conduct  a study,  design  and  fabricate  a 
cold  crucible  system  for  the  synthesis  and  crystal  growth  of  refractory 
semiconductors  and  oxides.  A theoretical  analysis  of  the  cold  crucible 
r«  lating  the  properties  of  the  material  to  the  electrical  characteristics 
the  cold  crucible  was  completed.  The  thermal  analysis,  which  includes 
tnermal  modeling  by  means  of  a computer  technique,  was  also  completed.  A 
detailed  analysis  of  the  RF  power  control  system  necessary  to  provide  the 
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desired  stability  for  crystal  growth  was  also  carried  out. 

A cold  crucible  assembly  suitable  for  the  growth  of  • i.vtory  oxide  crystal 
has  been  designed  and  constructed.  It  is  n s-mpl  , > -jggeo  modular  assenbly 
which  will  be  readily  adapted  in  the  future  ? > - i<‘C  atmosphere 

furnace  system. 

Using  a 50  Kw  (4-7  megahertz)  RF  power  supply  manufactured  bv  Fritz 
Hue t tinge r Elektronik  GmbH,  Freiburg,  German/,- li.i  '.old  crucible  assembly 
was  tested  and  kilogram  melts  of  zirconia  (melting  point  2690®C)  were 
successfully  contained. 

The  cold  crucible  assembly  developed  during  the  course  of  this  program  was 
delivered  to  Air  Force  Cambridge  Research  Laboraioi  ies  (LQP) . 
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‘"Theoretical  analysis  of  the  electrical  characteristics  of  the  cold 
crucible  system. 


It  was  found  that  a properly  designed  skull  melting  assembly 
has  no  effect  on  the  field  applied  to  the  charge.  The  skull 
merely  lowers  the  inductance  of  the  KF  coil  slightly  by 
excluding  the  field  from  the  volume  occupied  by  the  tube 
structure  and  adds  a small  resistive  loading  to  the  coil. 


The  power  delivered  to  the  molten  charge  by  a given  RF  field 
was  calculated  for  disc-shape!  charges  approximated  by  an 
ellipsoid  of  revolution.  These  calculations  were  carried  out 
for  both  large  and  small  field  penetration  skin  depths  and, 
for  an  ellipsoid  of  dimensions  8cm  x 8cm  x 2cm  in  a field  of 
1000  A/m,  the  maximum  power  required  is  several  kilowatts. 

In  considering  skull  melting  assemblies  of  varying  sizes,  our 
analysis  shows  that,  for  skin  depths  which  are  small  corap.'.red 
to  the  radius  of  the  skull,  the  power  required  is  generally 
proportional  to  the  square  of  the  radius.  However,  in  the 
case  where  skin  depth  is  large  compared  to  the  skull  radius, 
the  RF  power  must  be  transferred  at  a high  field  and  therefore 
lower  overall  efficiency.  In  general,  large  cold  crucible 
assemblies  appear  to  be  more  efficient  than  small  systems. 

< I 
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vstem. 


l:sing  computer  techniques  for  thermal  modeling  of  various  skull 
melting  configurations,  we  found  that  the  molten  region  tend 
to  take  the  shape  of  a flattened  body  of  revolution.  With  the 
source  of  heat  near  the  equatorial  plane  near  the  surface  <>* 
the  body,  heat  travels  inwards  and  in  lost  by  radiation  and 
conduction  along  the  way. 

if  a substantial  heat  loss  from  the  melt  surface  occurs,  the 
shape  of  the  melt  is  very  sensitive  to  the  amount  thus  lost. 
Therefore,  upper  radiation  shielding  is  vital  to  the  stability 
of  melt  tei*>erature. 

A detailed  analysis  of  the  RF  power  control  system  necessary  to  provide 
the  desired  stability  of  operation. 

A closed  loop  control  system  was  used  to  maintain  constant 
power  input  to  the  melt.  Initially,  RF  current  feedback 
control  was  explored  but  improved  control  of  RF  power  was 

achieved  hy  usir.-,  the  oscillator  plate  current  as  the  feedback 
signal. 

The  response  time  of  the  control  system  was  approximately  2 
seconds;  however,  the  themal  time  constant  of  the  molten 
charge  was  much  longer  than  this  with  the  result  that  melt 

temperature  fluctuations  due  to  external  disturbances  were 
relatively  small. 

'DesA‘';*  construction  of  a cold  crucible  assenfcly  suitable  for  the 
growth  of  refractory  oxides. 
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The  resultant  design  was  based  upon  tlu  guidelines  establ lslied 
by  the  theoretical  analysis  as  veil  as  first-hand  knowledge  ol 
the  cold  crucible  configurations  which  are  currently  operating 
with  considerable  success  in  the  U.S.S.R.  The  cold  crucible 
systea  developed  is  a rugged,  completely  metallic,  nodular 
assembly  which  will  be  readily  adaptable  in  the  future  to  a 
controlled  atmosphere  furnace  system;  indeed,  the  structure 
has  been  designed  to  withstand  external  operating  pressures 
up  to  100  atmospheres. 

"Testing  and  evaluation  of  the  cold  crucible  assembly  using  stabilized 
zirconia  (melting  point  2690°C). 


Over  thirty  (30)  melting  experiments  were  carried  out  in  the 
cold  crucible  assertly.  The  power  supply  used  was  a 30  kilowatt 
high  frequency  RF  generator  (manufactured  by  Fritz  Huettlnger 
Zlektrooik  GmbH,  Freiburg,  Germany)  tuned  to  operate  at 
approximately  4 megahertz. 


Using  kilogram  charges  of  zirconia  powder,  densification  and 
melting  could  be  achieved  with  applied  RF  power  levels  as  low 
as  6 kilowatts.  At  15  kilowatts  of  applied  power,  it  was 
possible  to  obtain  complete,  stable  melting  of  the  zirconia 
charge. 

Large,  semitransparent  columnar  grain  single  crystals  of 
stabilized  zirconia  (typically  3-4sm  on  edge  by  5-7ma  in 
length)  were  obtained  by  slow  cooling  of  melts  contained  in 
the  cold  crucible  assembly. 


“Development  of  the  criteria  for  integration  of  the  cold  crucible  assembly 
within  a controlled  atmosphere,  high  pressure  furnace  charter. 
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One  major  concern  was  the  effect  of  electrical  grounding  or. 
the  performance  of  the  cold  crucible.  It  was  found  that 
stable  melts  could  be  achieved  and  maintained  quite  readily 
whether  the  cold  crucible  was  grounded  or  not  to  the  RK  power 
supply. 

As  noted  above,  the  cold  crucible  assembly  has  been  designed 
for  operation  at  external  pressures  up  to  100  atmospheres. 

The  major  problem,  as  yet  unresolved,  involves  the  development 
of  a low-loss  coaxial  RF  feedthrough  designed  not  only  to 
provide  maximum  RF  transmission  at  megahertz-frequencies  but 
which  meets  the  requirements  for  safe  operation  at  high 
pressures  (approximately  100  atmospheres). 

"Upon  completion  of  the  program,  the  prototype  of  the  cold  crucible  assembly 
was  delivered  to  the  Air  Force  Cambridge  Research  Laboratories  (LQP) . 
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PREFACE 


Major  advances  in  the  development  of  cold  crucibles,  i.e.  water  coo Jed 
metal  crucibles  for  use  with  RF  induction  heating,  used  for  “skull 
melting",  have  taken  place  in  the  past  two  years  with  significant  impact 
on  refractory  materials  technology.  For  the  first  time,  it  is  possible 
to  contain  stable,  uncontaminated  melts  of  even  the  most  refractory 
materials  for  which  no  crucible  exists. 

Russian  scientists  at  the  Lebedev  Physical  Institute  in  Moscow  have 
developed  skull  melting  to  a point  where  Is  has  become  a reproducible 
production  operation.  Oxide  melts  at  temperatures  up  to  3000 °C  with 
weights  up  to  25  kilograms  are  being  used  in  the  production  of  a vide 
variety  of  strategic  materials  Including  mlque  laser  crystals, 
refractory  optical  elements  and  melt-cast  ceramics. 

The  new  process  Involves  direct  high  frequency  induction  heating  of  the 
material  contained  in  a water  cooled,  crucible-like  structure.  The  melt 
formed  is  contained  by  a sintered  shell  or  “skull"  of  Identical 
composition  so  that  the  problems  of  reaction  and  contamination, 
traditionally  the  most  severe  problem  faced  in  the  containment  of 
refractory  melts,  have  been  virtually  eliminated.  The  melt  can  then 
be  cast  or  recrystallized  using  standard  methods  of  crystal  growth. 
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J.O  INTRODUCTION 


The  principal  objective  of  the  program  was  the  design  and  development  **• 
of  a skull  melting  system  which  can  be  used  for  the  production  of  high 
purity  oxide  single  crystals.  Typically,  most  of  the  cold  crucibles  in 
operation  today  have  evolved  by  empirical  trial  and  error  procedures.  \ 
theoretical  analysis  of  the  cold  crucible  has  been  carried  out  to  relate 
the  properties  of  the  material(s)  to  be  melted  to  the  electrical  and 
thermal  characteristics  of  the  cold  crucible.  Based  upon  this  analysis 
we  have  established  basic  guidelines  for  the  design  and  construction  of 
cold  crucible  asseablies. 


One  of  the  major  problems  as  yet  unresolved  in  the  use  of  the  cold 
crucible  for  the  growth  of  single  crystals  is  related  to  the  stabilization 
and  precise  control  of  melt  temperature.  We  have  analyzed  the  critical 
control  parameters  of  the  cold  crucible  and,  based  upon  the  unique  RF 
power  control  techniques  which  have  been  developed  in  the  United  States, 
have  defined  a control  system  to  provide  the  high  degree  of  stability 
necessary  for  crystal  growing  applications. 

A cold  crucible  asseably  complete  with  control  system  suitable  for  oxide 
crystal  growing  operation  has  been  designed  and  constructed.  Testing  and 
evaluation  of  the  cold  crucible  asseably  and  control  system  have  been 
carried  out  to  explore  the  operating  techniques  suitable  for  crystal 
growing  applications.  Stabilized  zirconia  (melting  point  2690*C)  has 
been  used  to  demonstrate  the  very  high  temperature  operation  of  the  cold 
crucible  under  oxidizing  conditions.  Large  columnar  crystal  grains  of 
partially  stabilized  zirconia  have  been  recrystallized  from  stable  melts 
contained  in  the  cold  crucible. 

The  ultimate  goal  of  this  work  is  the  integration  of  the  cold  cr„:MV»e 
system  into  a crystal  growing  furnace  chamber  capable  of  operation  under 
eon t re led  atmospheres  at  pressures  up  to  100  atmospheres. 


I 


it  appears  that  no  difficulties  will  be  encountered  with  the  integration 
of  the  cold  crucible  inside  the  high  pressure  furnace  chaaber.  However, 
it  will  be  necessary  to  develop  a suitable  coaxial  high  frequency  RF 
transmission  line  for  the  furnace  which  will  satisfy  the  requireaents  of 
the  AS ME  Pressure  Vessel  Codes  as  well  as  the  requireaents  for  low  loss 
RF  power  transaission. 
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2.0  TECHNICAL  DISCUSSION 


2.1  INTRODUCTION 

During  the  course  of  this  program  we  have: 

1.  Theoretically  analyzed  the  electrical  characteristics  of  the  cold 
crucible  system.  The  problems  involved  and  the  conclusions  derived 
are  presented  in  Section  2.2  while  the  derivation  is  discussed  in 
Appendix  A. 

2.  Theoretically  analyzed  the  thermal  characteristics  of  the  cold 
crucible.  The  conclusions  are  presented  in  Section  2.3  and  the 
derivation  shown  in  Appendix  B. 

3.  Theoretically  analyzed  the  control  of  the  RF  power.  This  is 
discussed  in  Section  2.4  and  the  derivation  is  presented  in  Appendix 


Designed  and  constructed  the  cold  crucible  assembly.  Refer  to 
Section  2.5  for  review  of  the  design  criteria. 

Tested  the  cold  crucible  assembly  with  stabilized  zlrconia  (melting 
point  2690*C)  for  demonstration  purposes.  Operating  procedures  and 
results  arc  reviewed  in  Section  2.6. 

Investigated  the  criteria  for  integration  of  the  cold  crucible 
assembly  in  a high  pressure  furnace  chanber.  Section  2.7  reviews 
the  problems  involved. 


Delivered  a cold  crucible  assembly  to  the  Air  Force  Cambridge 
Research  Laboratories  (LQP). 
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2.2  THEORETICAL  ANALYSIS  OF  THE  ELECTRICAL  CHARACTERISTICS  OF  THKCOLD 
CRUCIBLE 

•tie  following  problems  involved  in  skull.  melting  of  refractory  materials 
have  been  theoretically  analyzed: 


‘The  effect  of  the  copper  cage  on  the  magnetic  field  applied  to  the 
refractory  charge. 

“The  resistive  cod  inductive  effects  of  the  cage  on  the  radio  frequency 
coil. 

°Ti e effect  of  the  skin  depth  of  field  penetration  on  the  total  power 
input  to  the  charge. 

°The  spatial  distribution  of  the  power  input. 

°The  resistive  loading  of  the  RF  coil  due  to  power  absorbed  by  the  charge. 

°The  RF  magnetic  field  required  to  melt  metal  particles  in  the  charge  to 
achieve  startup. 

°The  effects  of  scaling  the  size  of  the  cold  crucible  assembly  on  the 
applied  magnetic  field  as  related  to  skin  depth. 

It  is  found  that  a properly  designed  cage  has  no  effect  on  the  field 
applied  to  the  charge.  The  cage  merely  lowers  the  inductance  of  the  RF 
coil  slightly  by  excluding  the  field  from  the  volume  occupied  by  the 
rage  tubing  and  adds  a small  resistive  loading  to  the  coil. 


The  power  delivered  to  the  molten  charge  by  a given  applied  RF  field  is 
calculated  for  a disc-shaped  charge  approximated  by  an  ellipsoid  of 
revolution.  This  is  done  for  the  field  penetration  skin  depth  both 
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large  and  small  compared  with  the  smallest  semi-axis  of  the  ellipsoid. 
These  two  limiting  cases  converge  twoard  a common  maximum  value  o!  the 
power  when  the  skin  depth  in  each  case  approaches  the  size  of  the  semi- 
axis. For  an  ellipsoid  of  dimensions  8cm  x 8cm  x 2cm  in  a field  of  1000 
A/m,  this  maximum  power  is  several  kilowatts. 

The  calculations  show  that  the  heating  power  is  strongly  concentrated 
around  the  equator  of  the  ellipsoid,  especially  when  the  skin  depth  is 
•.mall  compared  with  the  thickness  of  the  ellipsoid.  An  expression  is 
derived  for  the  resistive  loading  of  the  RF  coil  by  the  induced  currents 
in  the  molten  charge. 

Startup  by  means  of  2mm  diameter  metal  spheres  embedded  in  the  powdered 
charge  is  found  to  be  theoretically  feasible  but  requires  a starting 
field  perhaps  10  times  larger  than  that  needed  to  maintain  the  melt. 

Our  analysis  shows  that  for  skin  depths  which  are  small  compared  to  the 
radius  of  the  skull,  the  power  required  is  generally  proportional  to  the 
square  of  the  radius.  However,  in  the  case  where  the  skin  depth  is  large 
compared  to  the  radius,  then  the  power  must  be  transferred  at  a higher 
fiiJLd  and  therefore  reduced  overall  efficiency.  The  conductivity  of 
Zr02  * is  approximately  10  Mho/cm  at  its  melting  point  and  approxi- 
mately .35  Mho/c*  at  about  1300®C  with  a reasonably  linear  relationship. 
At  a frequency  of  3.5  megacycles  the  skin  depth  is  0.85cm  at  its  melting 
point  and  5cm  at  1300°C.  Therefore  the  skull  constructed  with  a diameter 
of  approximately  7cm  will  couple  very  easily  at  the  melting  point  but 
will  exhibit  reduced  coupling  efficiency  at  1300°C.  If  the  material  to 

be  melted  has  a conductivity  vs.  temperature  profile  as  shown  in  Figure 
(3) 

1 for  AI2O3  , the  dramatic  change  in  conductivity  (over  7.  orders  of 
magnitude)  shows  that  control  of  thin  material  at  or  near  the  melting 
point  is  quite  critical.  Once  the  alumina  has  solidified,  the  coupling 
efficiency  will  be  so  poor  that  it  cannot  be  remelted. 
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2.  $ THERMAL  ANALYSIS  OF  THE  COLD  CRUCIBLE 


The  thermal  analysis  of  the  feed  powder  and  the  melt  when  heated  by  tin 
induced  RF  currents  and  cooled  by  the  walls  of  the  container  was  completed. 
Ue  find  that: 


°The  molten  region  tends  to  take  the  shape  of  a flattened  body  o'" 
revolution;  it  acts  as  a fin,  with  the  source  of  heat  in  a ring  near  the 
equatorial  plane  near  the  surface  of  the  body.  Heat  travels  inwards  and 
is  lost  by  radiation  and  conduction  along  the  way. 

'if  a substantial  heat  loss  from  the  top  surface  of  the  melt  occurs,  the 
shape  of  the  melt  is  very  sensitive  to  the  amount  thus  lost. 

°There  should  be  a large  amount  of  powder  to  insulate  the  melt  on  its 
underside,  although  the  exact  amount  is  not  critical.  This  means  that 
the  molten  zone  should  be  confined  to  the  upper  portion  of  the  cooled 
container. 

Thermal  models  have  been  developed  for  typical  skull  melting  systems  and 
the  effects  of  heat  sources  and  sinks  upon  melt  shape  have  been  profiled 
by  computer  techniques  and  are  discussed  in  Appendix  B. 


2.4  CONTROL  OF  RF  POWER 


Stable  control  of  the  RF  power  is  required  for  successful  crystal  growth. 
There  Is  a basic  tendency  toward  instability  when  induction  heating  is 
used  directly  to  melt  semiconductors  and  refractory  oxides.  Therefore 
a closed  loop  control  system  was  used  to  maintain  a constant  power  input 
to  the  melt.  Initially  RF  current  feedback  control  was  explored  but 
improved  control  of  RF  power  into  the  contained  melt  was  achieved  by  the 
use  of  the  oscillator  plate  current  as  the  feedback  signal. 


1 
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After  melting,  the  balance  between  heat  generation  and  heat  lost; 
r*ri»e«v;ses  results  In  a stable  system.  TIu?  responae  time  ol  the  voni  nd 
system  war;  about  2 seconds.  The  thermal  time  wins t an l of  the  rh.irge  was 
much  longer  than  this,  with  the  result  that  temperature  I J actual  Ions  <iu< 
to  external  disturbances  were  relatively  small. 

Using  a 50  kilowatt  output  high  frequency  RF  power  supply,  we  would 
obtain  fields  stronger  than  10^  amperes  per  meter,  which  were  sufficient 
to  melt  kilogram  charges  of  stabilized  zirconia  and  maintain  them  in  the. 
fully  molten  state. 


2.5  DESIGN  AND  CONSTRUCTION  OF  THE  COLD  CRUCIBLE  ASSEMBLY 


Based  upon  the  electrical  and  thermal  analyses  which  have  been  carried 
out,  we  completed  the  design  of  a cold  crucible  assembly  which  is  shown 
in  Figure  2.  The  copper  structure  has  been  designed  for  high  pressure 
operation  and  integration  with  the  ADL  crystal  growing  furnaces  presently 
in  use  at  the  Air  Force  Cambridge  Research  Laboratories.  The  cold 
crucible  assembly  is  shown  in  Figures  3 and  4. 


2.6  TESTING  OF  THE  COLD  CRUCIBLE  ASSEMBLY 


During  the  course  of  the  testing  and  evaluation  phase  of  the  program, 
over  thirty  (30)  melting  experiments  were  carried  out  in  the  cold 
crucible  assembly.  The  power  supply  used  was  a 50  kilowatt  high 
frequency  RF  power  supply  manufactured  by  Fritz  Huettinger  Elektronik 
r.mbK,  Freiburg,  Germany,  tuned  to  operate  at  approximately  4 megahertz. 
The  cold  crucible  in  operation  with  a 1 kilogram  charge  of  stabilized 
zirconia  is  shown  in  Figure  5. 

To  set  up  the  cold  crucible  for  operation,  the  assembly  is  first  wrapped 


-15- 


Llu>  I ulios).  The* 
Cot*  1 {*»><,  wiiter  is 
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minute. 


the  skull  was  filled  with  stabilized  zirconiu  powder  and  approximately 
10  grams  of  zirconium  metal  chips  (1/8"  diameter)  were  burled  in  the 
center  of  the  powder.  The  RF  power  (approximately  4 megahertz)  was  then 
turned  on  and  increased  to  an  indicated  plate  voltage  of  approximately 
4-5  kilovolts.  Within  10  minutes  the  zirconium  chips  would  become 
incandescent  and  heat  the  surrounding  powder  charge.  As  the  zirconia 
charge  heated,  the  plate  voltage  dropped  and  the  current  increased, 
indicating  more  efficient  coupling  due  to  the  increase  in  electrical 
conductivity  of  the  heated  zirconia.  As  the  RF  coupling  increased,  the 
powder  would  densify  and  melt  and  it  was  necessary  to  add  more  powder 

to  the  molten  charge. 

While  densification  and  melting  could  be  initiated  with  RF  power  levels 
as  low  as  6 kilowatts,  applied  RF  power  to  excess  of  15  kilowatts 
produced  better  results. 

Figure  6 shows  a fused  stabilized  zirconia  (Zircoa  Corporation)  charge 
as  removed  from  the  cold  crucible  assenbly.  Note  that  the  outer 
sintered  shell  of  the  charge  (approximately  lmm  thick)  has  retained  the 
shape  of  the  multi-tube  structure  of  the  cold  crucible.  The  top  of  the 
molten  charge  shows  the  effects  of  continuous  out-gassing  during  the 
course  of  the  fusion  operation.  No  external  radiation  shield  was  used 
to  reduce  the  severe  radiation  losses  from  the  upper  sui face  ot  the 
molten  charge,  and  a solid,  translucent  crust  formed  above  the  melt 
which  had  to  be  repeatedly  punctured  to  permit  the  addition  of  feed 

powder. 

Figure  7 shows  Che  initiation  of  welting  using  approximately  6 kilowatts 
of  applied  RF  power  to  the  charge.  Rote  the  sintered  dark  wass  near  the 
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hot  turn  where  the  zirconium  charges  were  placed.  It  is  also  evident  that 
once  heating  starts,  the  last  place  to  heat  is  the  center  region.  Figure 
6 shows  Sample  3 with  an  Input  power  of  7 kilowatts;  clearly  more 
sintering  and  melting  is  evident.  Figure  9 shows  Sample  7 with  an  input 
of  9 kilowatts  power.  First  evidence  of  melting  is  now  obviously  visible 
near  the  top  of  the  charge. 


It  was  evident  that  some  form  of  radiation  shielding  was  required  to 
reduce  thermal  losses  from  the  incandescent  melt  and  to  minimize  the 
hard  crust  formation  above  the  melt  surface.  A zlrconia  ceramic 
radiation  shield  was  positioned  on  top  of  the  cold  crucible  structure. 
However,  it  was  noted  that,  with  the  RF  load  coil  positioned  close  to 
the  top  of  the  cold  crucible,  the  hot  ceramic  radiation  shield  was 
heated  inductively,  thus  reducing  the  total  power  available  to  fuse 
the  charge  and  destroying  the  radiation  shield. 


Figure  10  is  an  example  of  fusion  (charge  M-2)  with  the  coil  near  the 
top  of  the  cold  crucible.  Note  that  while  better  melting  was  observed, 
crust  formation  was  still  appreciable.  In  this  run  the  radiation 
shield  was  destroyed.  In  Figure  11  (Sample  M-3) , the  coil  was 
substantially  lowered,  showing  significantly  less  crust  formation; 
the  shield  remained  intact 

The  decision  was  made  to  reduce  the  size  of  the  charge  to  approximately 
700  grams  (thus  lowering  the  melt  level)  and  to  lower  the  RF  coil  so 
that  the  lowest  turn  was  positioned  at  the  bottom  of  the  cold  crucible. 
Subsequent  experiments  using  ceramic  radiation  shielding  (Norton  sintered 
zirconia  bubble  plate)  proved  to  be  quite  satisfactory  and  melts  could 
be  readily  maintained  with  significantly  reduced  crust  formation. 

Figure  12  (Sample  M-9)  was  run  using  Y2O3  stabilized  zirconia  and  cooled 
to  room  temperature  over  1 hour.  The  sample  had  a distinct  yellow  color, 
but  large  columnar  transparent  grains  were  evident. 
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Figures  13  and  14  show  the  results  of  Run  M-10.  A 700  gram  zircon! a 
powder  charge  (Johnson  Matthey  Chemicals,  Ltd.,  99.99%  purity)  was  used 
with  10  wt  X yttria  (Rare  Earth  Products,  Ltd.,  99.991  purity)  added  to 
prevent  the  destructive  phase  change  (cubic  to  tetragonal  to  inonoclinio) 
which  occurs  during  the  high  temperature  recrystallization  of  zirconia. 

The  charge  was  melted  with  15  kilowatts  of  app3ied  RF  power  (estimated 
loss  via  radiation  is  less  than  2 kilowatts).  The  melt  was  cooled 
slowly  over  a period  of  1 hour  and,  upon  removing  the  recrystallized 
charge  from  the  cold  crucible,  large,  semitransparent  coluiwiar  crystals 
of  zirconia  were  revealed.  Upon  visual  examination  of  the  colutroar 
crystals  (typically  2- 3mm  on  edge  and  5-6tmn  in  length)  translucent 
regions  were  dispersed  throughout  the  transparent  grains  indicating  that 
the  high  temperature  cubic-phase  was  not  fully  stabilized;  it  appears 
that  more  than  10  wt  X yttria  should  be  added  to  the  zirconia  charge  to 
fully  stabilize  the  cubic  form. 


With  regard  to  the  use  of  automatic  RF  power  control,  it  was  found  that 
the  automatic  control  system  ceald  no i.  be  utilized  effectively  during 
the  heating  or  cooling  phases  ot  the  melting  operation.  At  the  critical 
temperature  where  the  skin  depth  of  the  charge  was  approximately  chat  of 
the  cold  crucible,  the  RF  coupling  efficiency  would  change  drastically, 
thus  causing  the  control  circuit  to  overdrive  the  RF  generator.  Once 
the  stable  melt  had  been  achieved  (using  manual  control  methods)  the  KF 
power  control  system  proved  to  be  very  effective  in  maintaining  a stable 
power  input  to  the  fused  mass.  At  the  termination  of  the  experiment, 
the  cooling  operation  was  carried  out  by  manually  reducing  the  input 
power. 


2_._7_ [ NTECRATION  OF  THE  COLD  CRUCIBLE  ASSEMBLY  WITH  A HIGH  PRESSURE  FURNACE 

To  date,  all  of  the  melting  experiments  with  the  cold  crucible  assembly 
have  been  carried  out  in  air  (as  shown  in  Figure  5)  and  there  are  obvious 
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advantages  to  using  this  assembly  within  a controlled  atmosphere,  high 
pressure  furnace  chamber. 

One  major  concern  was  the  effect  o f electrical  grounding  on  the  ncrl ormance 
of  the  cold  crucible.  The  cold  crucible  has  been  operated  both  with  and 
without  grounding  to  the  RF  power  supply  and  stable  melts  could  be 
achieved  quite  readily  in  either  case. 

The  cold  crucible  assembly  was  designed  initially  so  that  the  support 
stem  could  be  sealed  into  the  base  port  of  the  standard  ADL  Model  MP,  HP 
or  HPCZ  Furnaces.  An  0-ring  sealed  support  structure  (to  adapt  the  cold 
crucible  to  these  furnaces)  has  been  designed  to  permit  operation  at 
pressures  ranging  from  10  **  torr  to  100  atmospheres.  The  cold  crucible 
assembly  could  easily  be  mounted  in  a fixed  position  at  the  base  of  the 
furnace  to  facilitate  Czochralski  crystal  growing  experiments.  Alter- 
natively, the  stem  of  the  cold  crucible  assembly  could  be  attached  to  a 
lower  drive  mechanism,  thus  permitting  vertical  movement  of  the  assembly 
for  Bridgman-type  crystal  growing  operations. 


In  order  to  fuse  large  melts  (of  the  order  of  1 kilogram) , it  trill  be 
necessary  to  add  feed  powder  to  the  cold  crucible  assembly  during  the 
course  of  the  melt-down  operation.  A suitable  powder  storage  vessel 
and  feed  mechanism,  capable  of  operating  at  gas  pressures  up  to  100 
atmospheres,  have  been  designed  for  glass-melting  applications  prior  to 
the  initiation  of  this  program  and  should  prove  to  be  useful  in  resolving 
this  operating  problem. 

:'he  major  difficulty  we  can  foresee  in  integrating  the  cold  crucible 
assembly  within  a controlled  atmosphere  furnace  chamber  involves  the 
low- loss  transmission  of  high  frequency  (megahertz  range)  RF  power  into 
the  sealed  metal  chamber.  It  will  be  necessary  to  design  a suitable 
coaxial  RF  feedthrough  which  will  not  only  provide  maximum  RF  transmission 
efficiency  but  meet  the  requirements  for  safe  operation  at  high  pressures 
(approximately  100  atmospheres). 
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FIGURE  3 ADL  COLD  CRUCIBLE  ASSEMBLY  - SECTIONED 
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FIGURE  8 ZIRCONIA  CHARGE  FUSED  WITH  APPROXIMATELY  7 KW  RF  POWER 
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FIGURE  10  ZIRCONIA  CHARGE  FUSED  WITH  APPROXIMATELY  II  KW  RF  POWER 
(RF  COIL  POSITIONED  NEAR  THE  TOP  OF  THE  CHARGE) 


FIGURE  12  ZIRCONIA  CHARGE  FUSED  WITH  APPROXIMATELY  15  KW  RF  POWER 
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FIGURE  14  COLUMNAR  GRAINS  OF  ZIRCON I A OBTAINED  ON  SLOW  COOLING 


l.t)  W INCLUSIONS 

A simple,  rugged  cold  crucible  assembly  hau  been  developed  which  can  he 
used  lo  contain  stable  kilogram  size,  refractory  oxide  melts  at 
temperatures  up  to  3000#C  in  air. 

“Theoretical  analyses  of  the  thermal  and  electrical  characteristics  of 
cola  crucible  operation  have  been  carried  out  and  verified  experimentally. 
3ased  upon  these  analyses  we  have  established  the  basic  guidelines  for 
the  design  and  construction  of  cold  crucible  assemblies. 

cThe  critical  control  parameters  of  the  cold  crucible  have  been  developed 
and  we  have  defined  an  RF  power  control  system  necessary  to  provide  the 
high  degree  of  stability  for  crystal  growing  operations. 

°The  cold  crucible  assembly  has  been  designed  and  constructed  to  withstand 
external  pressures  up  to  100  atmospheres  and  we  have  established  the 
criteria  for  Integration  of  the  cold  crucible  assembly  within  a high 
pressure  furnace  chamber. 
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4.0  RECOMMENDATIONS 


Theoretical  analysis  of  cold  crucible  melting  and  processing  indicates 

that  minimal  problems  will  be  encountered  in  developing  large-scale  cold 

crucible  systems.  There  is  already  dramatic  evidence  of  the  potential 

o> 

industrial  applications  of  this  new  technology;  Russian  workers  have 

successfully  contained  refractory  oxide  melts  weighing  up  to  23  kilograms 

and  molten  glass  charges  up  to  250  kilograms  in  weight  are  now  being 

'6) 

nrocessed  in  France.' 


'lo  date  most  cold  crucible  operations  are  carried  out  in  air.  Under 
controlled  atmosphere  conditions,  this  new  process  may  well  be  the  answer 
to  the  time  honored  problems  of  working  with  ultra-refractory  materials 
such  as  ferrites,  semiconductors,  chalcogenide-glasses , metal  alloys  and 
cermets  whose  melts  are  sensitive  to  ambient  atmosphere  composition  or 
pressure. 

It  appears  that  few  difficulties  will  be  encountered  with  the  integration 
of  the  cold  crucible  system  inside  . high  pressure  furnace  chamber. 
However,  it  will  be  necessary  to  develop  a suitable  coaxial  high 
frequency  RF  transmission  line  for  the  furnace  which  will  satisfy  the 
stringent  requirements  of  high  pressure  operation  and  minimal  RF  power 
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APPENDIX  A 


THEORETICAL  ANALYSIS  OF  THE  ELECTRICAL  CHARACTERISTICS  OF  SKULL  MELTING 


A.l  EFFECT  OF  THE  CAGE 

We  asauae  chat  the  cage  consists  of  vertical  copper  tubes  parallel  to  the 
axis  of  the  coll  as  shown  In  cross  section  In  Figure  A.l.  The  effect  of 
the  bottom  of  the  cage  Is  small  and  will  be  neglected  here. 


CROSS  SECTION  OF  SKULL  MELTING  CAGE 
FIGURE  A.l 

If  the  coll  has  n turns  and  a length  9.  and  carries  an  alternating  current, 
1,  the  rms  current  per  unit  length  is 


nl 


(A/m) 


(1) 


and  the  axial  magnetic  field  is 


H 


fw 


(2) 
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At  s requoncies  In  the  range  1-10  megahertz,  the  skin  depth  oi  penetration 
-i  the  field  into  the  copper  tubes  of  the  cage  is  only  about  / r.  10  *i». 
.'he  copper  therefore  essentially  excludes  the  field  and  act*:  like  an 
almost  perfect  diamagnetic  material.  Under  these  condition:;  each  tube 
has  a circulating  current  per  unit  length  equal  to  -i  (A/m). 

The  energy  stored  in  the  field  is 


E-—  H 2 £ (A  - A ) (Joule) 
2 o coil  cage 


where  the  A*s  are  the  cross-sectional  areas  of  coil  and  cage  and  where 
= Ur.  x 10  ^ Henry/m  and  is  the  permeability  of  free  space.  From  (2) 


the  expression  becomes 


p n2l2  (A  - A ) 
r o coil  cage 

2 ft 


The  energy  can  also  be  expressed  In  terms  of  the  net  self- inductance,  L, 
of  the  assembly  by  the  relation 


* - 2U2 


On  equating  (4)  and  (5)  we  obtain  for  the  inductance 

(A  - A ) 

L = nzPo  — — -£ — (Henry) 


The  cage,  therefore,  reduces  the  inductance  of  the  coil  by  a factor 
^'coil  ^cage)/A 

The  cage  also  adds  some  series  resistance  to  the  coil.  The  circumferential 
current,  i,  in  a tube  of  the  cage  is  distribution  in  depth,  x,  into  the 
copper  according  to  the  approximate  formula 

J = (7) 

(l 


J 


where  J is  the  current  density  (A/m),  and  is  the  skin  depth,  given  l>y 
the  rr let  ion 


JUZ  ..  /"  2 


tt  u f a 

o 


VI  ‘DO 
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( H) 


Here  f is  the  frequency  (Hz)  and  a is  the  conductivity  of  copper  (mho;;/m) 


The  power  dissipated  per  unit  area  of  tube  surface  is  then 


W” 


(9) 


The  total  power  for  N tubes  of  height,  h,  and  diameter,  D,  is 


u - N TtD  h l2  n2  Nir  Dh  l2 
2o6  * 2o£26 


(10) 


ana  the  effective  series  resistance  introduced  into  the  coil  is  thetsfore 


n2  Nit  Dh 

R ■ — „T*r 


(U) 


We  have  assumed  that  h < £.  If  the  tubes  are  longer  than  the  coil,  h,  in 
(11) , is  replaced  by  SL. 


A. 2 HEAT  INPUT  TO  THE  CHARGE 


We  assume  for  mathematical  convenience  that  the  charge  has  the  shape  of 
•m  oblate  spheroid,  which  is  a good  approximation  to  the  ac’ual  disc-like 
■;nape.  The  equation  of  the  surface  is 


02) 


where  a and  c are  the  semi-axes  of  the  spheroid,  and  z and  r are  vertical 
and  horizontal  cylindrical  coordinates. 
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The  distribution  of  induced  current  in  the  spheroid  depends  strongly  on 
the  skin  depth,  6,  in  the  molten  material.  The  two  Halting  cases, 

6 » c and  6 « c,  are  mathematically  tractable,  while  the  intermediate 
case,  S c,  is  too  difficult  for  analytical  treatment.  The  limiting 
cases,  however,  serve  to  bracket  the  intermediate  case  and  give  considerable 
insight  into  the  general  problem. 


A. 2.1  Large  Skin  Depth 

When  6 » c,  the  uniform  magnetic  field,  Hq,  produced  by  the  coil, 
penetrates  the  charge  with  little  distortion.  Under  these  conditions 
the  induced  circumferential  electric  field,  E^,  is  given  by  the  single 
circuital  emf  equation 

2*t2Eq  * - “ (irr2  B)  - - itr2  uq  ju>  Hq  (13) 

where  ~r2  B is  the  instantaneous  magnetic  flux  through  a circle  of  radius, 
r,  and  w is  the  angular  frequency. 

The  associated  current  density  is 

J9  * ° Ee  = - 2 OUo  “ r Ho  (U) 

The  rate  of  heat  generation  per  unit  volume  at  the  radius,  r,  is 

- ilJ9l2  - £ °“c2  “2  r2  Ho7  (15> 

The  power  per  unit  area  projected  on  the  (x,y)  plane  is,  from  (12), 

W"  = 2z  W’"  = £ op  2 m2  H 2 c r2  /l  - (lb) 

It  is  to  be  noted  that  W"  is  zero  at  both  r = 0 and  r = a and  is  maximum 
at  r = a/ffi  * 0.816a. 
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Ih*-  total  heating  rate*  for  the  whole  ellipsoid  is 


W 


I W 2nrdr 
0 


2v 


op  ‘ or’  a1*  c tl  ■’ 
lio  <> 


< I /) 


From  (8),  this  can  also  be  written,  in  terms  of  instead  oi  a,  as 


W = 


p u a '*  c H 
4"  o o_ 

15  6 ■ 


(18) 


As  a numerical  example  we  take  the  following  values : 

p «=  4 x 10  ^ ti/m 

O z: 

f ■ 10  x 10  Mtiz 

a * 0.04m 

c * 0.01m 

6 = 2c  = 0.02m 

H = 1000  A/tn 

o 


Then  V * 4000  watts. 

The  value  chosen  for  6 corresponds  to  e = 63  mho/m  and  is  probably  near 
the  limit  of  validity  of  the  theory.  The  value  of  Ho  corresponds  to 
that  for  a 5-tum  coil  of  length  10cm  carrying  a current  of  2C  amperes. 


A_._2 . 2 Small  Skin  Depth 

The  case,  4 <<  c,  may  not  occur  in  practice  for  the  range  of  conductivities 
likely  to  be  found  for  molten  refractory  materials  without  the  use  of 
f regencies  well  in  excess  of  10  megahertz.  However,  it  is  still  of 
valu*  to  determine  the  power  distribution  in  the  charge  for  this  limiting 
ca-.e . 


Like  the  cage  tubing,  the  spheroid  acts  like  a nearly  perfect  diamagnetic 
body  when  <■  « c..  The  current  distribution  on  the  surface  of  the  spheroid 
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is  such  that  the  external  magnetic  field  Is  the  same  as  that  produced  by 
a uniformly  magnetized  body  of  the  same  shape.  The  intensity  oi 
magnetization,  a,  (dipole  noment  per  unit  volume)  is  related  to  the 
applied  field,  Hq,  of  the  coil  by  the  formula 


m = 


H 

o 

1-N 


09) 


where  N is  the  demagnetization  factor,  which  depends  on  the  ellipsoid 
shape.  For  a sphere  N = 1/3.  For  an  oblate  spheroid  with  semi-axes, 
a,  a,  c,  the  factor  is  given  by^ 

1 + v2  i _i 

N = ( - ^ ) (1  - ± tan  1 k)  (20) 

where 

k - ff-l  (21) 

Table  A.l  shows  how  N varies  with  c/a. 


TABLE  A.l 


DEMAGNETIZATION  FACTOR  N 


c/a 

N 

0.0 

1.000 

0.2 

0.750 

0.4 

0.588 

0.6 

0.478 

0.8 

0.3% 

1.0 

0.333 
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We  now  apply  the  boundary  condition  that  the  tangential  component  of  H is 

continuous.  Nov  the  field,  H , just  outside  the  ellipsoid  is  known  to  be 

tangential,  as  Indicated  in  Figure  A. 2.  Therefore,  just  below  the  surface, 

(H.)  » H , as  shown  in  the  vector  diagram  in  Figure  A. 3 where  the 

l tan  e 

vector  sua  of  H and  H equals  the  internal  field,  H. . 

e n i 


MAGNETIC  FIELD  VECTORS 
FIGURE  A. 3 


From  the  diagram  we  obtain  the  result 

H » H.  sin  1>  * - a sin  ♦ (25) 

e l 

where  the  angle,  *,  is  defined  in  Figure  A. 2. 


On  substitution  for  tn  from  (19)  we  find  that 


U 

e 


H sin  ♦ 
o 

1-N 


(26) 


The  surface  current,  i,  in  A/a  is  equal  to  the  tangential  magnetic  field. 
It  is  therefore  given  by 


i 


H sin  <P 
o _ 

1-N 


(27) 
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!o  evaluate*  sin  ♦ , we  represent  the  ellipsoid  by  the  parametric  equations 


r = a cos 


(2  8) 


7.  = c sin  <P 


(29) 


Then 


, dz 

sin  ^ = “— 

ds 


(30) 


where  ds  is  an  element  of  arc  on  the  allipsoid  at  the  point  of  tangency 
in  Figure  A. 2 and  is  given  by 


ds  “ /dr'’  + dz'* 


(31) 


From  (28)  and  (29) , 


dr  *»  - a sin  d>  d$ 


(32) 


dz  * c cos  y dd> 


(33) 


Therefore 


ds  = d$/a2  sin2  d>  + c2  cos2  i> 


( ;4) 


On  substituting  for  dz  and  ds  from  (33)  and  (34)  into  (30),  we  obtain 


sin  ^ = 


c cos  4> 


/a2  sin2  <J>  + c2  cos2  o 


(35) 


rherc-fore,  from  (27), 


i-N 


c cos  <p 


< !(»  ■; 


/a2  sin2  $ + c2  cos2 


The  total  power  is 
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W * R i 2 rds 
s 


\T 

i 


(37) 


(38) 


On  substituting  for  i,  r,  ds,  and  from  (27),  (28),  (34)  and  (38)  into 
(37),  we  find 


W = 


ac  H /2 
o 


cos  d 


(1-N) 


a ..  c H*) 

— * a(l-N)  " (W> 


- 2 a sin  + c cos 


where 


F(x) 


-~m  1o«. 


(40) 


(x  -i) 


Table  A. 2 shows  how  the  second  factor  in  (39)  depends  on  a/c.  It  can  be 
seen,  that  the  power  is  quite  constant  for  c/a  between  0.1  and  1.0  (for 
fixed  a).  For  very  flat  ellipsoids  the  power  increases  slowly. 


TABLE  A. 2 

GEOMETRICAL  FACTOR  IN  EQUATION  39 

c F(~) 
c 


c/a a(l-N) 


0.0001 

24,020 

0.001 

18.101 

0.01 

12.534 

0.1 

8.190 

0.2 

7.510 

0.25 

7.465 

0.4 

7.564 

0.6 

8.121 

0.8 

8.744 

1.0 

9.425 
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The  radial  distribution  of  the  power  is  also  of  interest.  From  (17)  we 
find,  on  counting  the  power  on  both  top  and  bottom  of  the  ellipsoid,  chat 


dW 

dr 


A 11  li  C 
s_  o 

* O-N) 


where 


<5> 


~7>  “111  - <"-)] 

d'  a 


(41) 


(42) 


Table  A.  3 shows  how  G varies  with  r/a  for  the  case  c/a  = 1/4.  The  last 
column  m the  table  gives  the  fraction  of  the  power  within  a radius  r. 

It  is  seen  that  the  power  is  strongly  concentrated  at  the  equator  of  the 
ellipsoid. 


TABLE  A. 3 

RADIAL  DISTRIBUTION  OF  POWER  FOR  c/a  = 1/A 


L 

H 

G<f> 

w(H 

W(a) 

0.0 

.cooo 

.0000 

0.1 

.0010 

0.2 

.0083 

.0001 

0.3 

.0296 

0.  A 

.0757 

.0021 

0.5 

.1650 

0.6 

. 3317 

.0121 

0.7 

.6532 

0.8 

1.3A92 

.05A9 

0.9 

3.A09A 

1.0 

1.0000 

For  c/a  = 1/4  the  total  power,  from  (39)  and  Table  A.  2,  is 


W 


7.A7 


H -V 
o 

of- 


7.A7  H V- 
o 


(A3) 
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*V  ••/ji<  rt  that  tli *s  exprcfiKion  will  remain  approximately  vaJld  loi 

i > 12  = 0.005m.  At  f - 10  megahertz,  this  correspond:.  t<>  >•  1000  inlio/m 

For  a - 0.04m  and  = 1000  A/m,  as  before,  we  find  from  (4  i)  that 

W = 2400  watts  which  is  to  be  compared  with  the  value  of  4000  watts 
calculated  on  the  6 » c assumption  for  4 = 0.02m. 

On  substituting  for  o from  (8)  we  find  that  (43)  takes  the  form 

W = 3.73  u u a2  H 2 6 (44) 

o o 

which  shows  that  W tends  to  zero  for  very  small  6,  for  constant  w.  The 
limiting  expressions,  (44)  and  (18),  therefore  show  that  as  6 varies 
from  small  to  large  values,  due  to  variation  of  the  conductivity,  W 
passes  through  a maximum  which,  for  the  numerical  values  chosen  above, 
probably  lies  between  2000  and  4000  watts  and  occurs  for  6 * c = 0.01m. 

T*-.  is  to  be  noted  that  for  constant  o,  W increased  with  increasing  w. 
as  can  be  seen  from  (43) . 

A. 3 COIL  RESISTANCE  DUE  TO  LOAD 


On  substituting  for  Hq  from  (2)  into  (44) , we  can  express  W in  terms  of 
the  current  1 in  the  coil  by  the  formula 


W 


4.78  p u a2  c n2  l2 
o 

*.?- 


(45) 


The  load  therefore  acts  like  a resistance 

4.78  p <*>  a2  6 n2 

\ ^ °__2 (6  « c) 


Likewise,  from  (18)  we  obtain 


R 


L 


St 


4,  “o  - a“  c n* 

15  E2!2 


(6  >>  c) 


(46) 


(47) 
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A.  4 STARTUP 

in  order  to  start  the  inductive  heating  one  must  usually  adu  f.orae  high 
conductivity  material  such  as  particles  o i the  same  metal  as  is  present 
in  the  composition  of  the  charge  itself-  Typically,  the  metals  used 
are  available  in  the  form  of  random-sized  small  chunks,  flake.,  or  chips 
and  for  purposes  of  our  calculations  we  have  assumed  the  particles  to  hd 
spherical  with  radius  a and  conductivity  a.  In  the  high  frequency  field, 
iio>  the  sphere  acts  as  a diamagnetic  body  with  a demagnetization  factor 
N = 1/3  (see  Table  A.l).  From  (19)  the  effective  intensity  of  the 
magnetization  of  the  sphere  is 


m 


K 


(48) 


The  total  magnetic  moment  of  the  sphere  is 


M = — ua^m  = 


- 2ira3  H 


(49) 


For  the  special  case  of  the  sphere,  unlike  that  of  a spheroid,  the 
moment,  M,  acts  as  though  it  were  located  at  the  center  of  the  sphere. 
The  tangential  magnetic  field  at  the  surface  of  the  sphere,  due  to  M, 

is 


<V 


M sin  9 1 „ , 

~T^s-  - - 7 H0  sln  0 


(50) 


dipole 

where  9 is  the  angle  between  the  radius  vector  and  the  z-ax.is. 


Since  the  applied  field,  H , gives  a tangential  component,  - H sin  >' , 

o o 

at  the  surface  of  the  sphere,  the  total  tangential  field  is 


H = - ^ H sin  0 

t 2 o 


'/>!) 


The  circumferential  current  per  unit  length  on  the  sphere  is  therefore 
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J 


i.  = - .,  II  «in  0 

$ / o 


since  it  equals  the  discontinuity  in  H at  the  surface  of  the  sphere. 

The  surface  resistance  Rg  is  given  by  (38).  On  multiplying  this  by 

i 2 and  Integrating  over  the  area  of  the  sphere,  we  obtain  for  the  power 
9 

input  to  the  sphere 

it  3 H 2a2 

W = / <2^5) <}H02  8i°2  e>  2lTa2  sin  eda  “ fg~  (53) 

It  is  of  interest  to  consider  that  this'  power  is  balanced  by  blackbody 
radiation  from  the  surface  of  the  sphere  at  absolute  temperature,  T , 
with  the  surroundings  at  temperature , T.  Then 

3irH  2 

4ira2  o (T4  - T 4)  * ~~~  (34) 

OO  Oo 


or  3 2 

•.«-  - V>  - -d-  (55) 

-8  -2  -4 

where  o * 5.672  x 10  W m deg  is  the  Stefan-Boltzmann  constant, 
o 

(55)  shows  that  the  equilibrium  temperature  is  Independent  of  the  radius 

of  the  sphere.  For  the  case  of  a field,  * 1000  A/m,  as  before,  and 

with  0 = 5.0  x 10~5  mho/m,  a (T4  - T 4)  = 1500  W/m2,  and  T « 431  K. 

o o 

The  metal  sphere  would  therefore  not  melt  for  the  field  assumed.  However, 
fo a ten  times  greater  field  of  10,000  A/m,  which  is  not  unreasonable, 
we  get 


0 (T4  - T 4)  « 150,000  W/m2 

o o 


T = 1276  K 


which  will  melt  certain  metals. 
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A more  appropriate  assumption  for  skull  melting  is  that  the  metal  sphere 
i;  embedded  in  a powder  of  thermal  conductivity  K.  Then  the  orjuil  ibriuai 
temperature  is  given  by 


T - T 


For  K - 0.1  watt/(m  deg  K) , a - 1 x 10  m,  = 10,000  A/m,  and  the 
other  constants  the  same  as  before,  ve  find  that  W = 1.88  watts  and 
T = 1600°K.  The  particle  would  therefore  raise  the  temperature  of  the 
surrounding  powder  to  a value  where  inductive  heating  of  the  powder 
would  occur. 


A. 5 SCALING  LAMS 


It  is  on  considerable  interest  to  determine  the  scaling  laws  for  magnetic 
field  Hq  and  the  power  input  W when  the  geometry  of  the  skull  and  the 
charge  are  changed  by  some  linear  factor.  To  find  the  scaling  laws  we 
assume  that  the  rate  of  heat  loss  from  the  molten  charge  to  the 
surroundings  is  given  approximately  by  the  usual  formula 


h A (T  - T ) 
mm  o 


where  h is  the  heat  transfer  coefficiant,  A is  the  surface  area  of  the 

n 

melt,  T^  is  the  temperature  of  the  melt,  and  is  the  temperature  of 
the  surroundings.  Since  h,  T^  and  Tq  are  constant  for  a given  material, 
(A?)  can  be  written  as 


where  a,  as  above,  is  the  major  semi-axis  of  the  spheroidal  melt,  it  is 
assumed  that  the  shape  of  the  melt  remains  fixed,  so  that  c/a  is  constant. 
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•nd»-r  equilibrium  conditions  the  total  power  W supplied  to  Jin-  molten 
cnarge  equals  the  rate  of  heat  loss. 


Therefore 


W = W. 


Loss 


('i9'» 


In  the  case  where  the  skin  depth  6 is  large  compared  with  a,  we  find  from 
(18),  (50)  and  (51)  for  constant  to  and  o. 


H 

o 


_1 

3/2 

a 


(60) 


This  means,  for  example,  that  if  the  geometry  is  scaled  down  by  a linear 

factor  of  2,  H (the  high  frequency  field)  has  to  be  increased  by  a 
o 

factor  of  2/2  = 3.83.  A linear  factor  of  4 requires  an  eightfold  increase 

in  H . 
o 


Jn  the  other  hand,  in  the  case  of  a highly  conductive  melt,  (31)  along 
with.  (50)  and  (51)  , shows  that 


H = constant  (61) 

o 

From  (50)  and  (51)  the  scaling  law  for  power  input  is  simple 

W * a2  (62) 


whatever  the  skin  depth  may  be.  However,  in  the  case  of  large  f-,  since 
the  power  must  be  transferred  at  a higher  value  of  H , and  therefore 
Higher  coil  current,  the  coil  losses  are  larger  in  this  case  and  the 
efficiency  is  therefore  lower. 
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A. 6 RELATION  TO  PREVIOUS  WORK 


The  foregoing  theoretical  analysis  and  the  following  thermal  analysis  and 
power  control  analysis,  together  constitute  a more  detailed  study  of  the 

(3) 

process  outlined  in  a recent  paper  by  Aleksandrov  et  al.  The  results 
of  the  study  generally  support  the  conclusions  of  the  Russian  autho.-s  and 
serve  to  sharpen  their  criteria  for  stable  operation  of  the  skull  melting 
technique. 


i 

■j 

S 

j 

-* 


i 

i 
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THEORY  OF  THE  THERMAL  CHARACTERISTICS  OF  SKULL  MELT  I NC 

It  has  been  established  already  chat  the  currents  produced  in  an  oblate 
spheroid  by  an  induction  heater  are  largely  concentrated  around  an 
equatorial  ring  which  is  located  near  the  surface.  While  the  current 
density  as  a function  of  position  has  been  derived  rigorously  only  for 
the  oblate  spheroid,  the  results  for  this  case  give  us  much  insight  and 
enable  us  to  anticipate  what  will  probably  happen  in  the  case  of  a more 
complicated  body  of  revolution  whose  radius  is  large  compared  with  its 
thickness  in  the  axial  direction  and  which  is  an  electrical  conductor. 

Qualitatively  speaking,  we  may  say  that  the  currents  In  the  body  will 


establish  themselves  in  such  a way  as  to  exclude  most  of  the  field 
impressed  by  the  induction  heater  coil.  The  value  of  B produced  by  a 
ring  of  current  depends  mainly  on  the  current,  and  not  so  much  on  the 
radius  of  the  ring  (in  a solenoid,  the  value  of  H is  independent  of  the 
radius  of  the  coil),  whereas  the  total  flux  pr>duced  is  the  product  of 
the  area  of  the  ring  and  the  average  value  of  B within  the  ring. 

Therefore,  in  the  case  of  a conductor  which  is  a flattened  convex  body 
of  revolution,  one  may  expect  that  the  induced  current  will  be  largely 
concentrated  in  a ring  of  the  largest  possible  radius  which  is  contained 
in  the  body.  It  is  probably  best  to  leave  this  kind  of  qualitative 
argument  in  its  present  rather  vague  state,  and  simple  to  conclude  that 
it  is  of  interest  to  study  the  shape  of  the  isotherm  at  the  melting 
temperature  in  a region  filled  with  a powder  of  low  thermal  conductivity, 
which  contains  a ring  source  and  which  is  fairly  closely  contained  by  a 
vessel  (represented  by  an  isotherm  at  300°K)  maintained  at  room  temperature. 


gp/astwMwaggta 


One  way  to  approach  the  problem  is  to  take  the  shape  of  the  300K  isotherm 
• rot:  the  drawing  of  the  Lubes  containing  the  coolant  and  to  .solve 
rigorously  the  boundary  value  problem  thus  defined.  However,  the  physical 
parameters  Involved  (for  example  the  thermal  conductivity  of  the  powder, 
which  will  depend  upon  its  packing  and  any  sintering  that  may  take  place) 
are  very  poorly  known,  so  it  is  not  justified  to  use  such  rigorous 
procedure.  Instead  we  may  employ  a simpler  technique  and  obtain  at  least 
a preliminary  understanding  of  the  effects  on  the  molten  zone  that  will 
be  produced  by  variations  in  the  main  features  of  the  system.  Without 
great  difficulty  some  important  results  can  thus  be  obtained  which  will 
serve  as  guide  posts  in  the  experimental  program. 


In  doing  the  simplified  problem  we  are  at  liberty  to  place  sources  and 
sinks  of  various  strengths  at  such  locations  as  we  choose,  provided  that 
these  locations  are  outside  the  region  of  interest,  which  we  will  take 
to  be  bounded  by  the  300K  isotherm.  Of  course  the  ring  source,  which 
is  used  to  model  the  induction  heating,  must  be  within  the  melt.  We  make 
use  of  external  sources  and  sinks  to  model  the  loss  of  heat  from  the  top 
and  to  control  the  shape  of  the  300K  isotherm. 


It  appears  that  the  300K  isotherm  must  pass  close  to  the  ring  source  in 
that  melt  in  order  to  give  the  Isotherm  at  the  melting  point  its  expected 
flat  shape.  Also,  we  want  the  300K  isotherm  to  be  close  to  the  melt  at 
the  top  (to  model  the  loss  of  heat  through  the  aperture  used  to  extract 
the  crystal)  and  otherwise  to  resemble,  at  least  roughly,  the  shape  of 
the  cup  formed  by  the  tubes  containing  the  coolant. 

In  Figure  ii.l  we  show  what  we  have  taken  to  be  the  "standard"  arrangement 
for  the  300K  and  2318K  isotherms  (the  melting  point  of  alumina  is  2318K) . 
"he  positions  of  the  sources  and  3 sinks  are  shown  and  the  drawing  is  to 
scale,  as  indicated.  The  arrangement  pictured  in  Figure  B.l  was  obtained 
after  a number  of  mathematical  experiments;  one  lesson  that  was  learned 
very  early  is  that  the  300K  sink  on  the  axis  has  a drastic  effect  on  the 
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shape  of  the  melt,  and  that  the  lower  .surface  oi  the  melt  needs  to  be 
well  insulated.  The  melt  is,  in  effect,  acting  as  a fin.  with  the  heat 
traveling  towards  the  center  and  being  lost  along  the  way.  The  insulation 
along  the  under  side,  which  serves  to  keep  the  heat  loss  to  a minimum, 
can  be  obtained  by  making  sure  that  the  molten  zone  is  in  the  top  part  of 
the  cooled  container,  so  that  a thick  blanket  of  powder  is  left  unmelted 
beneath  the  2318K  isotherm. 

We  will  now  illustrate  the  sensitivity  of  the  shape  of  the  molten  region 
to  changes  in  some  of  the  system  parameters  by  means  of  figures  which 
show  the  two  standard  isotherms  and  two  slightly  distorted  isotherms. 

The  distorted  isotherms  are  obtained  by  solving  the  problem  with  some 
changes  in  the  boundary  conditions. 

In  Figure  B.2  we  show  two  pairs  of  isotherms  to  demonstrate  the  sensitivity 
of  the  shape  of  the  melt  to  a variation  in  the  position  of  the  upper  part 
of  the  300K  isotherm  from  the  standard  position.  The  lower  part  of  the 
300K  isotherm  has  been  held  nearly  in  place.  It  is  clearly  seen  that  the 
top  surface  of  the  melt  moves  further  than  the  300K  isotherm,  while  even 
the  bottom  part  suffers  a significant  displacement.  It  is  of  interest  o 
note  that  the  changes  seen  are  due  mainly  to  a change  of  50  watts  in  the 
strength  of  the  sink  on  the  axis  (from  300W  to  250W) . 

In  Figure  B.3  we  show  two  pairs  of  closely  spaced  isotherms.  The  lower 
part  of  the  300K  isotherm  has  been  caused  to  move  about  3n»>  below  its 
standard  position  (on  the  axis),  but  the  corresponding  variation  in  the 
position  of  the  2318K  isotherm  is  only  about  1mm  at  the  bottom  and  nearly 
nothing  at  the  top  of  the  melt.  Thus  it  is  seen  that  the  shape  of  the 
molten  zone  is  not  very  sensitive  to  variations  in  the  geometry  of  the 
bottom  of  the  cup,  provided  that  the  blanket  oi  insulation  is  present. 

We  also  determined  the  sensitivity  of  the  shape  of  the  molten  zone  to  a 
change  in  the  input  power  (supplied  by  the  induction  heater)  by  raising 
the  strength  of  the  ring  source  from  3150W  to  3300W,  while  the  shape  of 
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the  300K.  isotherm  was  kept  essentially  the  same  as  the  standard.  The 
result  is  shown  in  Figure  B.4,  which  indicates  that  tlu*  shape  of  the  melt 
is  only  slightly  affected  by  the  change  In  power  input.  This  is  probably 
due  to  the  fact  that  the  distance  between  the  side  of  the  molten  zone  and 
the  cold  container  is  very  small,  so  a minute  change  in  this  dimension  is 
all  that  is  needed  to  cause  most  of  the  added  power  to  flow  directly  into 
the  container  walls  near  the.  plane  z = 0. 

Under  certain  circumstances  it  might  be  desirable:  to  let  the  molten  mass 
cool  slowly,  of  its  own  accord,  without  penetrating  the  layer  of  powder 
that  lies  over  the  melt.  In  such  a case  on<  might  obtain  a symmetrical 
melt,  as  shown  in  Figure  B.5. 

So  far,  we  have  discussed  only  the  case  of  molten  alumina.  Another 
material  of  interest  is  zirconia,  whose  melting  point  is  298SK.  In  Figure 
B.6  we  show  a set  of  sources  and  sinks  and  the  corresponding  isotherms  at 
2938K  and  300K.  Because  of  its  higher  temperature  the  upper  surface  of 
the  pool  of  molten  zirconia  may  be  expected  to  lose  much  more  heat  by 
radiation  than  is  the  case  with  alumina.  To  allow  for  this,  the  strength 
of  the  point  sink  has  been  set  at  800W,  and  the  strengths  of  the  other 
sources  and  sinks  adjusted  in  an  attempt  to  produce  a suitable  configura- 
tion of  the  2983K  and  3G0K  isotherms.  While  the  amount  of  mathematical 
experimentation  with  the  model  of  the  higher  melting  point  material  has 
been  very  limited,  we  expect  that  it  would  be  necessary  to  separate  the 
ring  source  slightly  further  from  the  main  ring  sink  in  order  to  produce 
a 2588K.  isotherm  that  is  similar  in  appearance  to  th»2  2318K  isotherm  of 
Figure  B.l.  From  an  operational  point  of  view,  this  means  that  the 
frequency  of  the  induction  loeater  must  be  set  so  as  to  obtain  a largo 
skin  depth  in  the  melt. 

We  expect  that  the  sensitivity  of  die  shape  of  the  zirconia  melt  to 
changes  in  the  external  configuration  will  be  qualitatively  the  .•ame  as 
in  the  case  of  the  alumina  melt  already  discussed. 
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APPENDIX  C 


TEMPERATURE  DISTRIBUTION  DUE  TO  A RING  SOURCE 


In  this  section  ve  derive  the  formula  for  the  temperature  (in  steady  state) 
at  a point  in  an  infinite  medium  which  contains  a source  of  heat  in  the 
form  of  a ring.  Ve  assume  the  thermal  conductivity  of  the  medium  to  have 
the  fixed  value  K (watts  cm  deg  C ) , so  the  temperature  nust  satisfy 
Laplace’s  equation 


72T  = 0 


and  furthermore  we  require  that  the  temperature  be  zero  at  infinity. 


In  accordance  with  these  assumptions  we  find  the  temperature,  dT,  at  a 
field  point,  P,  whose  coordinates  are  (x,y,z)  produced  by  a point  source 
of  strength  dQ  watts  located  at  P*  with  coordinates  (x1 ,yf,z*)  to  be 


dT 


_iSL 

4*Ks 


(1) 


where 

s = Ax-x’)2  + (y-y')2  + (z-z’)2  (2) 

By  adding  up  the  contributions  due  to  many  such  sources  we  can,  at  least 
in  principle,  find  the  temperature  due  to  any  distribution  of  heat  sources. 

In  Figure  C.l  we  show  a ring  source  of  radius,  a,  located  in  the  z-plane 
and  a field  point  at  an  arbitrary  point.  It  is  evident  from  the  diagram 
that 

x’  » a cos  0 

y'  * a sin  0 

z’  = 0 (3) 
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Now  suppose  Che  source  to  dissipate  Q*  watts  cm  uuiloriuiy  along  Its 
circumference.  Then 


dQ  = Q*  ad8  (4) 

since  ad  is  the  eleaent  of  arc.  Therefore  we  find  from  the  above  equati 
equation 

dT  - -~~~~de --; (5) 

4irK/(x  - a cos  6) 2 + (y  - a sin  9) 2 + z2 

or,  by  integration  of  (5)  around  the  loop. 


T 


Q*a  ** d9 

0 J (x  - a cos  8) i + (y  - a sin  0)2  + z2 


(6) 


On  account  of  the  symmetry  of  the  problem  we  can  replace  x in  (6)  by  r 
and  replace  y by  0.  Thus 


T 


Sl±y  

0 /r2  + a2  + z2  - 2ar  cos  8 


(7) 


The  integral  in  the  right-hand  side  of  (7)  is  not  "elementary",  but  can 
be  evaluated  simply  in  terns  of  the  well  known  complete  elliptic  integral, 
K(m) , which  is  defined  by  the  formula 


K(m) 


0 /l  - m sin?  <J> 


First  make  the  substitution 
in  (7)  to  obtain 


6 


2$  + * 


T 


XS 


9l± 

2*K 


*/2 

I 


d$ 


-it/2  /(a  + r)2  + z2  - 4ar  sin2  $ 


(8) 


-9) 


Next  remove  the  factor  (a  + r)2  + z2  from  the  expression  whose  square 
root  is  indicated  and  note,  since  sin2  $ is  even,  that  the  range  of 
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APPENDIX  D 


CONTROL  OF  RADIO  FREQUENCY  POWER 


The  typical  current  feedback  control  system,  shown  in  Figure  D.l,  was 
used  to  control  the  output  of  the  radio  frequency  generator.  The  Leeds  & 
N'orthrup  H-AZAR  Recorder  was  used  as  a preamplifier,  the  input  summing 
point,  and  also  to  record  the  error  signal.  The  characteristics  time 
constant  (t^)  of  the  H-AZAR  is  0.5  second.  The  time  constants  of  the  RF 
generator  and  the  associated  saturable  reactor  are  approximately  this 
»ar-e  value  or  a little  longer. 

To  adjust  the  controller  to  the  process  the  integral  (reset)  action  (kc) 
of  the  Leeds  & Northrup  Series  80  CAT  Controller  was  made  as  large  as 
possible  in  order  to  obtain  accurate  steady  state  control,  while  the  gain 


(k^)  and  the  derivative  (rate)  action  (k^)  of  the  controller  were  adjusted 
to  obtain  stable  control  with  minimum  response  time.  The  resultant 
settings  were:  Reset  - 50  repeats /minute.  Proportional  band  = 35%,  and 

Rate  = .02  minute.  The  frequency  response  of  the  system  is  shown  in 
Figure  D.2.  The  open  loop  response  was  computed  from  the  measured  closed 
loop  response  values.  The  1/f  slope  in  the  low  frequency  end  of  the  open 
loop  response  curve  is  due  to  the  integral  action  of  the  controller.  Ihe 
65°  phase  margin,  close  to  the  minimum  value  considered  acceptable  in 
good  control  practice,  shows  that  the  adjustment  is  close  to  optimum. 


excessive  pickup.  Some  effects  of  RF  noise  were  still  noticeable  even 
after  our  best  efforts  to  Improve  the  system  grounding. 

Another  factor  which  interfered  with  closed  loop  control  oi  the  RK  power 
during  startup  was  the  dramatic  change  in  conductivity  of  zirconia  upon 
heating.  The  conductivity  rises  rapidly  with  temperature  and  there  is  a 
correspondingly  rapid  increase  in  power  absorbed  by  the  load.  At  the 
same  time  there  is  little  change  in  the  RF  current  because  the  melt 
occupies  only  a small  portion  of  the  volune  inside  the  heating  coil  and 
also  because  the  penetration  depth  is  still  large  compared  to  the 
dimensions  of  the  melt.  The  result  is  a tendency  toward  thermal  runaway. 

To  solve  this  problem  the  Huettinger  RF  generator  was  modified  to  allow 
the  oscillator  plate  current  to  be  used  as  the  feedback  signal.  The 
plate  current  increased  almost  in  direct  proportion  to  the  real  power 
load  on  the  RF  generator  and  is  a much  better  measure  of  the  power  into 
the  melt  than  RF  current.  Also,  this  method  of  control  was  not  sensitive 
to  RF  noise  pickup.  When  the  control  system  was  adjusted  for  plate 
current  feedback,  the  controller  settings  were:  Reset  * 100  repeats/ 

minute.  Proportional  band  * 100Z,  and  Rate  = .02  minute.  The  overall 
response  was  not  measured  but  seemed  to  be  about  the  same  as  with  the  RF 
current  feedback. 

The  last  factor  in  the  system  dynamics  is  the  thermal  transfer  function 
of  the  charge.  This  function  is  difficult  to  measure  but  we  can  estimate 
the  thennal  time  constant  (t^)  by  assuming  that  the  heat  loss  is  entirely 
due  to  conduction  and  using  the  relationship 

t3  - «hf  <» 

The  heat  capacity,  = 4.18  Ms  (Joules/degree)  in  which  M and  s are  the 
mass  and  specific  heat  of  the  charge,  respectively.  Typical  values  are 


-68- 


astossaBBSsgsiiassagaagBMaM £ J frr ssa  rr>r  v 


M = 700  gm 

s * 0.1  cal/gm°K 


| 


i 


% 


k 


K 


l 


which  result  in  C.  = 293  Joule/°K. 
h 

The  thermal  resistance,  dT/dW,  is  approximately 

AT 

dT/dW  « — <degrees/watt) 

in  which  AT  is  the  temperature  rise  of  the  melt  and  W is  the  total  input 
power. 

Evaluating  this  derivative  for  AT  = 2600°C  and  W = 13,000  watts,  we 
obtain  dT/dW  « .2  degrees  per  watt  per  degree  and  = 60  seconds.  Since 
this  is  much  longer  than  the  response  time  of  the  closed  loop,  temperature 
fluctuations  will  be  considerably  smaller  than  the  power  fluctuations  of 
the  system. 

One  final  consideration  is  the  maximum  available  power.  Power  is 
determined  by  the  magnetic  field  and  therefore  by  the  RF  current.  The 
current  is  determined  by  the  inductive  reactance  of  the  coil  and  the 
output  voltage  of  the  RF  generator.  Since 

»*r*  I = fe  L*“2  (2) 

decreasing  the  number  of  turns  in  the  coil,  n,  increases  the  current,  I, 

2 

as  n and  increases  the  magnetic  field,  H,  and  consequently  the  output 
power. 


Our  final  coil  design  is  for  a 4-tum  coil,  2"  long  and  2.1V*  mean  radius. 
The  coil  inductance  is  given  by 


L 


n2  r2 
10  + 9r 


42  x 2. IS2 

10  x 2 + 9 x 2.15 


1.9  microhenries 


**3) 
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The  t actor  lor  loss  of  inductance  due  to  the  cross-sectional  area  of  the 
cage,  (.Area  - Area  ; /Area  coil,  is  estimated  to  be  about  .8,  so 
chat  the  net  inductance  is  about  1.5  microhenries.  Assuming  that  the 
maximum  voltage  is  5,000  volts  and  at  a frequency  os  3.6  megahertz,  we 


I * — -_-r  - 150  amperes 

2n  x 3.6  x 10b  x 1.5  x 10 


4 x 180 

2"  x 2.54  x 10" 


= 11,600  amperes/meter. 


Substituting  this  value  into  (17)  from  Appendix  A, 


W * — a p A u.,?-  a**c  H 2 
15  o o 


with  the  other  parameters 


we  obtain 


* 50  mho/meter 

* 4 x 10  7 henry/meter 

= 19  x 10^  radians/second 

= .032  meter 

= .01  meter 


W = 24  kilowatts 


which  is  half  the  rated  capacity  of  the  RF  generator. 


In  practice  we  measured  20  kilowatts  plate  power  at  approximately  these 
conditions.  The  power  into  the  cooling  water  was  about  13  kilowatts. 

This  was  close  to  the  maximum  power  obtainable  in  a practical  sense  since 
at  higher  voltages  there  was  a tendency  toward  arcing  across  the  work 
coil. 
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